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Abstract 

We illustrate a method for improving Renormalisation Group improved pertur- 
bation theory by calculating the infra-red central charge of a perturbed conformal 
field theory. The additional input is a dispersion relation that exploits analyticity 
of the energy-momentum tensor correlator. 

1 Introduction 

The purpose of this letter is to describe a method for improving a Renormalisation 
Group improved perturbative calculation. We will study a two dimensional quantum 
field theory off criticality. In the two scaling limits, the infrared (IR) and the ultra- 
violet (UV), it becomes a conformal field theory characterised by the Virasoro central 
charge. The field theory we look at has the minimal models A4(m) as scaling limits. 
Our approximation method obtains a value for the infrared central charge Cjr using 
perturbative information around the ultraviolet limit of the theory. In |1|, [J, [| cm was 
calculated in perturbation theory in the limit when m —>■ oo. This is valid because 
then the UV and IR fixed-points are arbitrarily close in the coupling constant space. 
The renormalisation group (RG) eigenvalue y = ^-j- of the perturbation vanishes as 
m — > oo and the perturbation thus becomes marginal, and hence the field theory does 
not move away from the UV fixed-point. As well as for large values of m our ap- 
proximation also applies for small values, (m > 10), and we can test it against the 
exact result cjr = CM(m-i) when cjjv = CMim)- in the limit m — > oo we obtain the 
perturbative result reported in ^ f|. To improve upon standard perturbative tech- 
niques requires some additional input. For this we will exploit the analyticity of the 
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energy-momentum correlator (TT) in a complex scale parameter using it to construct 
a dispersion relation. In section 2 we describe this dispersion relation. Section 3 gives 
the results from perturbative CFT and the RG improvements, and in section 4 we 
illustrate our approximation with numerical results. 



2 The dispersion relation 



We will now construct a dispersion relation that relates the infra-red and ultra-violet 
behaviour of the energy-momentum tensor two-point function. The correlator of two 
energy-momentum tensors can be written using the Kallen-Lehman spectral represen- 
tation H as 

71 r a ~t f d2 P e^ pZ+pz) _ 4 _ F(\z\ 2 ) 
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where we use the usual complex variables z, z, and c(fi)dfi is the spectral density of 
the QFT which represents the density in degrees of freedom at the mass //. Integrating 
over p gives 

1 r°° - 

F(s) = — c//ic(/i)/i V ^((/xV /2 + 24/iVi)^o(^) + (8/i 2 s+48)X 1 (^)) (1) 
4o Jo 

where s G M+ and Kq, K\ are the modified Bessel functions. In H it was shown using 
the spectral representation that cyy = J °° dfj,c(fj,) and Cm = lim e ^ fo df/,c(jj,). From 
the properties of the modified Bessel functions it then follows from (|l|) that F(s) has 
the limits 

F(s) - ( Cuv I 01 S ° + ' (2) 
[ c ir tor s — ► oo. 

In we prove in detail that the expression (|I|) provides an analytic continuation from 
real positive values of s to the complex plane cut along the negative real axis, as might 
be expected from the analyticity of the Bessel functions themselves. From now on 
we take s to be in the complex plane with a thin wedge about the negative real axis 
removed (so — 7r + e < arg(s) < 7r — e). From ([l]) it also follows that F(s) has the limits 
(0) for all s — >• and |s| — ^ oo from the cut complex plane, so lini| s |__ >0O F(s) = cjr and 
F(0) = c uv . 

The following contour integral in the cut complex plane therefore vanishes by ana- 
lyticity of F(s) 

— [ ds—F(s) = — ( [ ds+ [ ds+ [ ds^]—F(s) = 0, (3) 
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where the contour C is given in figure 1. The contribution from the large circle C2, 
where s = r IR e ld for 9 G [n — e, — n + e], picks out the value (—1 + -) lini| s |_ >o0 F(s) = 
(— 1 + ^)cm in the limit rm 00. This is seen by writing the integral as an angular 
integral over 9, the limit rm — > 00 can then be taken before the integral as it follows 
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Figure 1: Integration contour C in the cut complex plane. 

from ([I]) that the integrand is bounded by a constant in the cut plane. For e small the 
infrared central charge can thus be calculated from an integral over the small circle Co 
and an integral over the contour C\ which we denote Cut(p) 

Cjk — 7t~t [ ds —F{s) + Cut(p) = I(p) + Cut(p). (4) 

* m J c s 

The exponential factor in the integrand suppresses the contribution of the cut for large 
positive values of p and lim^oo Cut(p) = so that cj R = lim^oo I(p). To show this we 
divide the \s\ interval into [r^y, 1] and and note that F(s) is finite in the unit 

disk hence the contribution in the lower region is bounded by the function k— for some 
k e M+. In the other region, using the properties of the Bessel functions, the integral 
is seen to be bounded by the expression k J dpc(p)p 7 / 2 J^ IR d re -P cos ( e )A> 3 / 4 e~^v / ^ sii W 2 
which is uniformly convergent in p, the details are shown in ||. 

For small s, F(s) and hence I(p) are determined by the ultra-violet behaviour. 
Close to the UV fixed-point we may describe this by perturbation theory with running 
coupling constant g(s), vanishing as s — > 0. We denote the perturbative approximation 
to the integral in (f|) by I n (p) where F(s) is approximated to n-th order by the pertur- 
bative expression F n (s). If analyticity is not spoilt by the perturbative expansion (as 
may be seen by inspection of the results in the next section) then simply substituting 
/„ for I in this limit would yield cjr ~ lim^oo I n (p) = lim^oo F n (s), which is just 
the perturbative estimate of the central charge, c* IR . Since the scale dependence of the 
theory can be absorbed in the scaling of the running coupling constant g(s) we can 
write F n (s) = $ n (g(s)), so that if g* IR is the first non trivial zero of the perturbative 
/^-function i.e. g(s) — > g} R for s — > oo, then lim p ^ 00 / n (p) = & n (gJ R ). However, we 
wish to improve on this result. Both I n (0) and 1(0) equal cuv- F° r small enough p, 
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I n {p) provides a good approximation to I(p), since the power series expansion of I(p) 
is controlled by the small s expansion of F(s) for which perturbation theory is good. 
For larger values of p, I(p) and I n (p) part company, tending to Cm and c* IR respec- 
tively. If Cm < c* IR and if the region where I n (p) is a good approximation to I(p) is 
large enough, then I n (0) — I n (p) will have a maximum before approaching its limiting 
value of cuv ~ c* IR . It is this maximum that we will use to provide a better estimate 
of cuv ~ c iRi since this occurs at the largest value of p for which I n (p) is a reasonable 
approximation to I(p) and the true value of cuv ~ c m is given by 1(0) — I(oo). As we 
will see this is the case for the minimal models considered in the next section. 



3 The perturbative calculation 

In the CFT M.(m) the primary field 0(i,3) is a relevant field which RG trajectory is a 
geodesic in the coupling constant space ||. The QFT 



S = S M ( m ) - A J d x 0(1,3) (a)- (5) 

therefore interpolates between the UV fixed-point (Ao = 0) given by A4(m) and the IR 
fixed-point which is given by the CFT M.(m — 1). This can be seen from perturbative 
arguments in the m — > oo limit |l|, ||, || and the general case is argued for in using 
thermo-dynamical Bethe ansatz methods. 

When Ao is small we can use perturbation theory to calculate quantities in the 
theory (Q). Calling the bare field 0(i,3) = 4> one gets to order A 



(0(x)0(O)) 



W^y) V + y 11 + ( °' 



where A(y) = r(i-y/2)^r(i+!/) = 1 + 0(y 3 ) and b(y) is the operator product expansion 
coefficient Cf, which can be calculated in the Coulomb gas representation of minimal 
models using the formulas in ||. b(y) is given by 

h f..\2 _ 16 (1-y) 4 f r(l+y/2) \ 4 f r(l-y/4) \ 3 ( T(l-y) \ 2 / r(l+3y/4) \ _ 16 , n (..\ 

U W) 3 (l-y/2)2(l-3y/4)2 \T(l-y/2)) \T(l+y/A) ) \T(l+y) J \T{l-3y/4) J S^^y). 

y = is the RG eigenvalue for the perturbation 0(i,3). With the renormalisation 
conditions ((p(x, g)(f>(0, g))\\ x \ =fM -i = p 4 , the renormalised correlator and /5-function, 

dfi 



f3(g) = p^j^-, becomes to the same order || 



p 4 / , 4ttv%)%)o 



(<Kx,g)<K0,g)) = + 

/%) = -y9-^b(y)g 2 A(y) + O(0 !i ). 
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Here g and (j>(x, g) are the renormalised coupling and field. The RG fixed-points, the 
zeros of the /^-function, are thus g uv = 0, g* IR = wA [ y ^ y) and therefore g E ( ^jfej , 0) 
as the theory (|5|) lies between the two scaling limits. The trace of the energy-momentum 
operator is the infinitesimal generator for scale transformations and it is given by [1| 
Q(x) = 2nf3(()(x). Then it follows from the Ward identities that 

d Ml dg 2 (T(z!, zi)T(z2, z 2 )) = ^d Zl d Z2 (Q{z 1 , z±)Q(z 2 , z 2 )) 

\ (6) 
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13(g) d Zl d Z2 {(j)(z 1 ,z 1 )(f>(z 2 ,z 2 )). 



Writing F(R) = 2z A (T(z, z)T(0, 0)) in terms of the dimensionless variable R = fi 2 zz, 
then (H) becomes 



dR 2 ~ K "' V " dR 2 

The limiting values of F(R) are given by @, using these the solution to (0) becomes 

fVf^ ir 2 g 2 R? fy(2-y)(3-y) n , f(2-y)(3-y) 

F(R) = c uv + -^r— ~ + 2tx A(y)b(y)g 1 



y - i V i - y 

5g (3y-4)(3y-6) 

3(|y - 1) 

The function F(s) is determined as F(s) = F(R)\p_ 1 s , g y The theory is thus fixed at 
the renormalisation point and all scale dependence is moved into the running coupling 
constant || 

y 

g(s) = 



wA(y)b(y)g ^ _ ^ 



We thereby get our 1 loop RG improved approximation of F(s) 

jpf n k 2 _ 2 f y (2-y)(3-y) 

Fits) = c uv + tt9 (« 



2 v ' \ y-l 

{2-y)(3-y) , (3y - 4) (3y - 6) 



+27rA(y)fo(y)£(s 



1-2/ 3(f y — 1) 



(9) 



To apply our approximation based on the dispersion relation of the previous section to 
calculate Ac = cuv ~ c m we now have to find the maximum of 



im-h{p)=c uv -^-. I ds—F^s). (10) 

Zm -'Co s 
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4 Results 



We first note that i*i(s) has the correct value in the perturbative limit y — > 0. From 
([5p it follows that Ac peri = cuv — c* IR = lim s ^ (cuv — -^i( s )) is given by 

a _ it 2 2 fy(2-y)(3-y) 

^Cp er t — ^ ly/Kj 



and the exact value is 

Ac ea: a rf = c(m) - c(m - 1) = — — =-= r = % r = — + 0(y 4 ). (12) 

K J K J m (m 2 -l) 2(2 -y)(4-y) 16 vy ; v ; 

To calculate Ac for finite m we have to calculate h(p). The running coupling @ can 
be rewritten as 

-/ \ 9^ g* IR \g\s y/2 _ 5 

= = TT^' 9 = i^f R ( -°°' 0) - (13) 

A variable change s > = s \g\ 2 /y in h (p) then leads to 

= h L ds ir F ^ = 2^? / c , = A&o I") 

where now all dependence of the renormalised coupling g is moved into p' = p\g\ 2 ^ y and 
the radius of C' : r' uv = r uv \g\ 2 l y as now g(s') = 1 ^yj/ 2 ■ The integral /i(0) — /i(p') 
can be calculated numerically and in figure 2 we have plotted it against log p' for the 
case of m — 14. For large p' it tends to the perturbative value 0.0025, but has a 
maximum value of 0.00367 which provides a better approximation to the true value of 
0.00440. The dashed line indicates how we expect 1(0) — I(p') to behave. We have 
approximated Ac for different values of m by the maximum of h(0) — I\(p') and we 
denote this approximation as Ac approx . The numbers Ac ap prox have been obtained by 
a numerical integration using a NAG mark 18 Fortran Library quadrature routine. In 
figure 3 the error in Ac approx and Ac pert compared with Ac exact is plotted against m. 
The errors are scaled with m(m 2 — 1) so that all points are distinguishable on the same 
plot. The figure shows that the error in the approximation Ac appTOX is more than a 
factor 2 smaller than the error in the perturbative value Ac per t in the region plotted. 

For values of m smaller than 11 perturbation theory begins to breakdown indicated 
by the RG improved result turning negative (violating unitarity). Our approximation 
whilst still positive for m close to 10 also becomes poorer since it is based on the RG 
expression. Also, note that both Ac approx and Ac per t approaches Ac exac t faster than 
the asymptotic value 3y 3 /16 in the limit m — > oo. 
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Figure 2: The numerical-result I\(0) — h(p') against logp', also plotted is the exact 
value Ac exact and the RG improved perturbative result Ac pert all for m = 14. The 
dashed line is the expected behaviour of 1(0) — I(p'). 
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Figure 3: {Ac exact - Ac approx )m{m 2 - 1) and (Ac exact - Ac pert )m(m 2 - 1) against to. 
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5 Conclusion 



We have described an approximation method to obtain the IR central charge of the 
minimal models valid for m > 10. The method improves upon the RG improved 
calculation by exploiting the analyticity of the correlator of the energy-momentum 
tensor. This analyticity is a property both of the exact correlator and the RG improved 
perturbative estimate of it, ||, which we have used in our calculation. 

Our approximation is correct in the perturbative limit m —>■ oo. For smaller m 
we obtained the values shown in figure 3 together with the RG improved perturbative 
result. This figure demonstrates that the approximation is significantly better than the 
RG improved perturbative result, e.g. for m = 18 it has a 7.8% relative deviation from 
the exact result whereas the RG perturbative result deviates by 18.0%. The analyticity 
in the complex scale parameter s of the energy-momentum tensor two-point function is 
in fact an ubiquitous property of correlation functions in quantum field theory, having 
its origin in the hermiticity of the Hamiltonian. Consequently we expect our approach 
to be applicable beyond the specific calculation we have used to illustrate it here, to 
other RG improved perturbative calculations of Green's functions. 
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